Geometric and dynamic perspectives on phase-coherent and noncoherent chaos 
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Statistically distinguishing between phase-coherent and noncoherent chaotic dynamics from time series is 
a contemporary problem in nonlinear sciences. In this work, we propose different measures based on recur- 
rence properties of recorded trajectories, which characterize the underlying systems from both geometric and 
dynamic viewpoints. The potentials of the individual measures for discriminating phase-coherent and nonco- 
herent chaotic oscillations are discussed. A detailed numerical analysis is performed for the chaotic Rossler 
system, which displays both types of chaos as one control parameter is varied, and the Mackey-Glass system as 
an example of a time-delay system with noncoherent chaos. Our results demonstrate that especially geometric 
measures from recurrence network analysis are well suited for tracing transitions between spiral- and screw-type 
chaos, a common route from phase-coherent to noncoherent chaos also found in other nonlinear oscillators. A 
detailed explanation of the observed behavior in terms of attractor geometry is given. 



PACS numbers: 

Oscillatory processes can be frequently observed in nat- 
ural and technological systems. Often, the corresponding 
dynamics is not strictly periodic, but shows more complex 
temporal variabiUty patterns characterized by a fast di- 
vergence of trajectories with arbitrarily close initial condi- 
tions [1-3 |. There are numerous examples of such chaotic 
oscillators for which long-term predictions of amplitudes 
and phases are not possible. Therefore, studying their 
phase dynamics has recently attracted particular inter- 
est, e.g., regarding the process of phase synchronization 
between different coupled systems |4, 5|. However, most 
existing methods suitable for this purpose require the ex- 
pUcit definition of an appropriate phase variable, which 
can become a non-trivial problem in the case of nonco- 
herent chaotic oscillations. Therefore, studying the phase 
coherence properties of chaotic systems has become an 
important problem in both theoretical and experimental 
studies |6|. In this work, we propose some methods based 
on the concept of recurrences in phase space, which allow 
studying complementary aspects of chaotic oscillators re- 
lating to the geometric structure of, and the dynamics on 
the attractor. Specifically, we derive a detailed character- 
ization of changes of the geometric structure of complex 
systems in phase space with varying control parameter, 
which accompany transitions from phase-coherent to non- 
coherent dynamics. 



I. INTRODUCTION 

In the last decades, the complexity of chaotic oscillators 
has been widely characterized by a variety of different quanti- 
ties inspired from nonlinear dynamical systems theory QH). 
Lyapunov exponents ||9l [TOll describe the characteristic time- 



scale associated with the finite-time exponential divergence 
of nearby chaotic orbits and, thus, relate directly to the pre- 
dictability horizon of the dynamics. Fractal dimensions and 
entropies measure the structural complexity of the underlying 
attractor, often based on concepts from information theory. 

In contrast to the aforementioned concepts, in many situa- 
tions, one is interested in explicitly characterizing the phase 
dynamics of the recorded nonlinear oscillations. However, 
depending on the structural properties of the chaotic oscilla- 
tions under study, it may be difficult to assign a well-defined 
phase variable to the observed dynamics. This problem pre- 
dominantly occurs in the presence of noisy oscillations, how- 
ever, also in the fully deterministic case, one frequently ob- 
serves oscillations without a distinct center of rotation in 
phase space, e.g., in the funnel regime (see Fig. [TJi) of the 
Rossler system ifTTl 

i = -y- z, 

y = x + ay, (1) 
i = 0.4 + z(a;- 8.5). 

In case of such noncoherent oscillations, the appropriate def- 
inition and analysis of the phase dynamics is challenging. 
Therefore, given the rising number of examples of real-world 
chaotic oscillators, the problem of automatically distinguish- 
ing between phase-coherent (PC) and noncoherent (NPC) 
chaos is of practical relevance. Traditionally, this problem has 
been considered by studying the phase diffusion properties of 
the system under study 1 5 1 . However, in order to apply this 
conceptual idea, an appropriate phase variable has to be de- 
fined in advance. 

In this work, we propose an alternative approach based on 
the recurrence properties of a dynamical system's trajectory 
in phase space for quantitatively characterizing whether or 
not an observed chaotic dynamics is phase-coherent. In con- 
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FIG. 1 : Two-dimensional projection of a part of the trajectory of the 
Rossler system (Eq. l[T}) in the (A) PC (a = 0.165) and (B) NPC 
(funnel) regime (a — 0.265). 



trast to the explicit study of phase diffusion, the corresponding 
concepts do not rely on an explicit definition of a phase vari- 
able. We emphasize that this fact has already motivated us- 
ing recurrence-based properties for studying synchronization 
processes of coupled NPC oscillators lfT2l [T3 1 and time-delay 
systems I14HT51. 

Generally, recurrence properties can be conveniently ana- 
lyzed by using recurrence plots (RPs) 1131 . originally intro- 
duced in the seminal work by Eckmann et al. which 
provide an intuitive visualization of the underlying temporal 
structures. For this purpose, one defines the recurrence matrix 
Rij as a binary representation of whether or not pairs of ob- 
served state vectors on the same trajectory are mutually close 
in phase space. Given two state vectors and (where 
i and j are time indices), this proximity is most commonly 
characterized by comparing the length of the difference vec- 
tor between Xj and Xj with a prescribed maximum distance e, 
i.e., 



^ l|x»-xj||), 



(2) 



where Q{-) is the Heaviside function and || • || a norm (e.g., Eu- 
clidean, Manhattan, or maximum norm). In this work, we will 
specifically chose the maximum norm for defining distances 
in phase space, since it has lower computational demands than 
other norms. However, the choice of a different norm would 
not change the presented results qualitatively. The proper- 
ties of RPs have been intensively studied for different kinds 
of dynamics |13|, including periodic, quasiperiodic |17 -19|, 
chaotic, and stochastic dynamics ||201|2TI . 

It has been shown that, among other features, the length 
distributions of diagonal and vertical structures in RPs can be 
used for defining a variety of measures of complexity, which 
characterize properties such as the degree of determinism or 
laminarity of the system Il22lj25l . The resulting toolbox of 
recurrence quantification analysis (RQA) has been widely ap- 
plied for studying phenomena from various scientific disci- 
plines lfT3l l26l . In this work, however, we will utilize some 
complementary conceptual approaches also based on RPs, 
which do not belong to the set of classical RQA measures. 
Based on Eq. (j2]i, we will discuss properties based on the re- 
cuiTence time statistics and so-called e-recurrence networks 
(RNs). The underlying methodological concepts are briefly 
described in Sec. |ll] and subsequently applied to two realiza- 
tions of the Rossler system in PC and NPC regime, respec- 



tively. Following the results obtained for this example, po- 
tential new statistical indicators for phase coherence based on 
the recuiTence properties of the underlying system are intro- 



duced in Sec. Ill and compared to other established as well 



as novel measures based on phase diffusion and Poincare re- 
turn times, respectively. Application to a complete bifurcation 



sequence of the Rossler system in Sec. IV demonstrates the 
feasibility of the recurrence-based approaches. The geometric 
consequences of the transition from PC (spiral-type) to NPC 
(screw-type) chaos and their impact on the recurrence prop- 
erties are discussed. As a second example, the behavior of 
the recuiTence based measures is illustrated for the Mackey- 
Glass system |27| in a parameter range including transitions 
between periodic and NPC chaotic behavior | 



II. METHODS 
A. Recurrence time statistics 

Complementary to RQA, another natural way for charac- 
terizing the recurrence properties of dynamical systems in 
phase space is statistically evaluating the distribution of re- 
cuiTence times (RTs), which has been applied to both chaotic 
and stochastic systems |29-32|. In contrast to return times 
with respect to a fixed Poincare surface, recurrence times re- 
fer to the time intervals after which the trajectory enters the 
e-neighborhood of a previously visited point in phase space. 
Gao et al. 1331 demonstrated that similar to some line-based 
RQA measures, characteristics based on the RT distributions 
p{t) can be used for detecting subtle dynamical transitions, 
which motivated using a corresponding approach for testing 
against stationarity L34ll35l . Besides their immediate impor- 
tance for studies on extreme events ll36l . recurrence times 
have also proven their potential for the estimation of dynami- 
cal invariants such as the information dimension |30| and the 
Kolmogorov-Sinai entropy |37|. 

Given a RP, RTs can be identified as the lengths of non- 
interrupted vertical (or horizontal, since the recurrence matrix 
is symmetric) "white lines" that do not contain any recurrence 
(i.e., no pair of mutually close state vectors). More precisely, 
such a white line of length r starts at the position in the 
RP if II38I 



1 if TO = —1, 

for TO e {0, , 

1 if TO = T. 



1}, 



(3) 



In order to see this, for all times k = j — 1, . . . , j + t, the 
values Xfc on the trajectory are compared with x^. Then, the 
structure given by Eq. ([3]) can be interpreted as follows: At 
time /c = j — 1, the trajectory falls into an e-neighborhood of 
Xi. Then, for fc — j, . . . , j + t — 1, it moves further away 
from Xi than a distance e, until at k = j + t, it returns to the 
£-neighborhood of x^ again. Hence, given a uniform sampling 
of the trajectory in the time domain, the length of the line is 
proportional to the time that the trajectory needs to return e- 
close to Xi. Going beyond the concept of first- return times, 
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the ensemble of all recurrences to the e-neighborhood of 
induces a RT distribution for this specific point. Combining 
this information for all available points in a given time se- 
ries (i.e., considering the lengths of all white lines in the RP), 
one obtains the RT distribution p(t) associated with the ob- 
served (sampled) trajectory in phase space. Hence, the length 
distribution p{l) of "white" vertical lines I in the RP not con- 
taining any recurrent pair of observed state vectors provides 
an empirical estimate of the distribution of RTs on the con- 
sidered orbit, which contains important information about the 
dynamics of the system under investigation. 



B. Recurrence network analysis 

Recently, different approaches have been proposed for 
studying basic properties of time series from a complex net- 
work perspective Il39l - l44l . Many existing methods for trans- 
forming time series into network representations have in com- 
mon that they define the connectivity of a complex network 
- similar to the spatio-temporal case - by the mutual proxim- 
ity of different parts (e.g., individual states, state vectors, or 
cycles) of a single trajectory ll44l |45]| . Among other related 
approaches, e-recurrence networks (RNs) and their quanti- 
tative analysis have been found to allow identifying transi- 
tions between different types of dynamics in a very precise 
way |i43]|Sli6l|42l|50l. 

In order to construct the RN, we re-interpret the recurrence 
matrix Rij, the main diagonal of which is removed for conve- 
nience, as the adjacency matrix Aij of an undirected complex 
network associated with the recorded trajectory, i.e., 

= - (4) 

where 6i,j is the Kronecker delta. The vertices of this net- 
work are given by the individual sampled state vectors on 
the trajectory, whereas the connectivity is established accord- 
ing to their mutual closeness in phase space. This definition 
of a complex network provides a generic way for analyzing 
phase space properties of chaotic attractors in terms of net- 
work topology B31 l50l . However, since the network topol- 
ogy is invariant under permutations of vertices, the statistical 
properties of RNs do not capture the dynamics on the attrac- 
tor, but its geometric structure based on an appropriate sam- 
pling. In this respect, we emphasize that since a single finite- 
time trajectory does not necessarily represent the typical long- 
term behavior of the underlying system, the resulting network 
properties depend - among others - on the length N of the 
considered time series (i.e., the network size), the probability 
distribution of the data, embedding [51 1, sampling |44, 52], 
etc. We choose the threshold e in such a way that the resulting 
RN has a fixed edge density (recurrence rate) of RR = 0.03 
unless otherwise stated explicitly. 

Although they primarily describe geometric aspects, the 
topological features of RNs are closely related to invariant 
properties of the underlying dynamical system P3l |45] l50l 
l53l . In model systems (e.g., Rossler and Lorenz systems), 
both local and global network properties have already been 



studied in great detail 1041 145] |50l |5T]. Among others, two 
particularly interesting local measures are 

(i) the local clustering coefficient Cy, which quantifies the 
relative amount of closed triangles centered at a given 
vertex v (i.e., at the associated point x^, in phase space) 
and gives important information about the geometric 
structure of the attractor within the e-neighborhood of 
V in phase space 1501 , and 

(ii) betweenness centrality by, which quantifies the fraction 
of all shortest paths in a network that include a given 
vertex v ||54| . In a RN, vertices with high 6„ corre- 
spond to regions with low phase space density that are 
located between higher density regions. Hence, by yields 
information about the local fragmentation of the attrac- 
tor ||45ll5Tl . Since in a complex network, the values of 
by may span several orders of magnitude, in the follow- 
ing we will consider log by as a characteristic measure 
for network topology. 

In a RN, both Cy and by are sensitive to the presence of un- 
stable periodic orbits (UPOs), but resolve complementary as- 
pects |51|. Specifically, in a continuous system, it is well- 
established that if a chaotic trajectory enters the neighbor- 
hood of an UPO, it stays within this neighborhood for a cer- 
tain time [55|. As a consequence, states accumulate along 
this UPO instead of homogeneously filling the phase space in 
the corresponding neighborhood (in particular if we consider 
UPOs of lower period), which results in a locally reduced ef- 
fective dimension that can be quantitatively characterized by 
Cy and measures derived from this quantity ISOl . 

In addition to the aforementioned vertex characteristics, 
several global network measures have already proven to dis- 
tinguish between qualitatively different types of behavior in 
both discrete and continuous-time systems ll43l l47r .49 1 . Ex- 
tending these previous results to different appearances of 
chaotic dynamics, we will consider four particular mea- 
sures 1 5 6. -.5 8 J as potential candidates for discriminatory statis- 
tics: 

(i) the global clustering coefficient C ll59l . which gives the 
arithmetic mean of the local clustering coefficient C„ 
taken over all vertices v, 

(ii) network transitivity T ll60ll6T1l . which is closely related 
to C (but gives less weight to poorly connected ver- 
tices 1 50 1) and globally characterizes the linkage rela- 
tionships among triples of vertices in a complex network 
(i.e., the probability of a third edge within a set of three 
vertices given that the two other edges are already known 
to exist) ||89l , 

(iii) the average path length C, which quantifies the average 
geodesic (graph) distance between all pairs of vertices, 
and 

(iv) the assortativity coefficient TZ |62|, which characterizes 
the similarity of the connectivity at both ends of all edges 
in the network (i.e., the correlation coefficient between 
the degrees of all pairs of connected vertices). 
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Network transitivity and average path length have already 
proven to provide an excellent discrimination between com- 
plex periodic and chaotic orbits in a two-parameter bifurca- 
tion scenario of the Rossler system [47 1. An analytical theory 
for computing the value of T from a known invariant density 
p{x) revealed a strong relationship to a certain concept of gen- 
eralized fractal dimensions [50 1. In this respect, high values 
of T indicate the presence of a lower-dimensional structure 
in phase space corresponding to a more regular dynamics. In 
contrast, the average path length behaves differently for dis- 
crete and continuous-time dynamical systems Il43ll45]l47l : for 
maps, more regular dynamics is characterized by low values 
of £, whereas the opposite applies to chaotic oscillators. 



C. Recurrence properties of phase-coherent and noncoherent 
Rossler systems 

As a simple continuous-time deterministic dynamical sys- 
tem that exhibits both PC and NPC chaotic dynamics, we 
first study the behavior of the RP-based concepts described 
in Sec. |II A| and |IIB| for the Rossler system (Eq. ([T}). In the 
following, we will use numerical simulations of this system 
for various parameters a, obtained with a fourth-order Runge- 
Kutta integrator with step width h = 0.01. The resulting tra- 
jectories have been down-sampled to = 10, 000 data points 
with a sampling interval of At = 0.2, which avoids strong ef- 
fects of trivial temporal correlations. 

In order to illustrate qualitative differences in the behav- 
ior of the RP-based indicators for PC and NPC dynamics, we 
consider the individual cases a ~ 0.165 (PC) and a = 0.265 
(NPC), respectively. A part of the resulting trajectories (pro- 
jected onto the (a;, y) -plane) is shown in Fig. [T] One clearly 
recognizes that the oscillations of the system have a well- 
defined center in the PC case, but no unique center for NPC 
chaos. 

The RT distributions obtained for both examples are qual- 
itatively different (see Fig. |2]i. Specifically, in the PC regime 
the lengths of time intervals without any recurrences are 
peaked around multiples of the basic period of oscillations, 
with a maximum at three full periods of the system [^8 1 (note 
the logarithmic units in Fig. |2]i. This indicates that in this 
regime, one distinct time-scale dominates the dynamics of the 
system. In contrast, in the NPC case, the distribution be- 
comes much more irregular, which indicates that a multiplic- 
ity of time-scales is relevant in the observed chaotic dynamics. 
However, since the complex structures in the RT distributions 
have not yet been expUcitly studied in previous work, it is not 
a priori clear which kind of statistical property (e.g., mean re- 
currence time or the corresponding standard deviation) can be 
used for distinguishing between both cases. Specifically, since 
the recurrence time is related to the mean period of oscilla- 
tions, its mean value varies considerably within the different 
dynamical regimes as the parameter a is changed. 

In contrast to the RT statistics, the local RN properties 
characterize higher-order features of the attractor geometry in 
phase space rather than dynamical aspects [45 [. While global 
network properties have been recently applied for automati- 
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FIG. 2: RT distribution p(r) with r — lAt (zoom for short times) 
for one realization of the Rossler system with (A) PC and (B) NPC 
chaos. The threshold e has been chosen to yield a recurrence rate 
RR = 0.03. 




FIG. 3: Color-coded representations of local RN properties ((A,B) 
local clustering coefficient C^, (C,D) logarithm of betweenness cen- 
trality log &„) for the Rossler system with (A,C) PC and (B,D) NPC 
chaos {RR — 0.03). In (C,D), black circles indicate vertices in 
poorly populated regions of phase space with 6„ < 1. 



cally discriminating between chaos and periodic dynamics in 
a complex two-parameter bifurcation scenario of the Rossler 
system f^7\, we suggest that local properties are able to char- 
acterize even more subtle structural changes of the system. 
For the two considered test cases. Fig. [3] shows the pattern of 
the local clustering coefficient Cv and betweenness centrality 
by in phase space. It is clearly visible that both measures char- 
acterize different aspects of attractor geometry [51 1, which re- 
sults in a correlation coefficient that is still significant, but not 
very large (Fig.|4|. Specifically, both measures are somewhat 
sensitive to the presence of UPOs, which are densely embed- 
ded in the chaotic attractor. However, while the corresponding 
direct relationship has been theoretically established only for 
Cy SO far in terms of an effective local dimension of the attrac- 
tor [50], by is no direct indicator for UPOs. 

Studying the full probability distributions of both local net- 
work measures in some more detail (Fig.|5]l, we observe clear 
differences between PC and NPC dynamics. Specifically, all 
distributions are at least bimodal (which is partially related to 
the presence of UPOs leading to locally increased clustering 
coefficients), whereas the bimodality is more expressed in the 
phase-coherent case. Together with the general finding that 
the maxima of the respective distributions do not differ con- 
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FIG. 4: Scatter plot between the RN measures Cv and log bv for the 
Rossler system with (A) PC and (B) NPC chaos {RR = 0.03). ps 
gives the values of the rank-order correlation coefficient (Spearman's 
Rho) between both quantities. 
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FIG. 5: Probability density function of the RN measures (A,B) Cv 
and (C,D) logfei, for the Rossler system with (A,C) PC and (B,D) 
NPC chaos. The different symbols represent the results obtained 
for the same trajectory with different choices of the recurrence rate 
(RR = 0.02 (□), 0.03 (*), and 0.04 (•)). 



siderably, this result motivates considering simple statistical 
properties of the distributions of C„ and log by for deriving 
novel indices for phase coherence. We will come back to this 
idea in the next section. 



III. QUANTIFYING PHASE COHERENCE OF CHAOTIC 
OSCILLATORS 

A. Phase and frequency of chaotic oscillators 

In order to numerically study phase coherence of chaotic 
oscillators, a reasonable definition of a phase variable is 
usually required first. While the derivation of optimum 
phase variables has been recently attracted considerable in- 
terest ||63,- 66J, we restrict our attention in this work to the 
standard analytical signal approach. Here, a scalar signal x{t) 
is extended to the complex plane using the Hilbert transform 



y{t) - -V.V. / \ ^ ' ds, 



(5) 



where V.V. denotes Cauchy's principal value of the integral, 
which yields the phase 



(j)(t) — arctan 



x{ty 



(6) 



We emphasize that this definition is appropriate for oscil- 
lations with a well-defined center in the origin of the (a;, y)- 
plane. Specifically, for PC dynamics, it is possible to find 
simple (linear) transformations of x and y (e.g., subtracting 
the mean) so that the oscillations are centered around the 
origin. In contrast, NPC dynamics is characterized by the 
non-existence of such a unique central point in phase space 
(cf. Fig.[Tj3). As a result, defining the phase in the above way 
leads to a variable that does not monotonously increase with 
time. Within the framework of phase synchronization anal- 
ysis, an alternative phase definition has therefore been pro- 
posed based on the local curvature properties of the analytical 
signal |[67H69l , i.e.. 



(j){t) = arctan 



dy{t)/dt 
dx{t)/dt' 



(7) 



We note that the proper evaluation of the derivatives in the 
latter equation may pose substantial numerical challenges, es- 
pecially in the case of (noisy) experimental data. 

The instantaneous frequency of a chaotic oscillator is de- 
fined as the derivative of the phase variable with respect to 
time. Averaging this property over time yields the mean fre- 
quency 



2it\ dt I 



(8) 



Since in the standard Hilbert transform-based defini- 
tion, the phase variable does not necessarily increase 
monotonously in time, we quantify this monotonicity in order 
to obtain a simple heuristic order parameter for phase coher- 
ence, which we will refer to as the coherence index 



CI 



lim — 

T->oo T 



e(-<^(i)) dt 



(9) 



with (\){t) = d4>{t)/dt. 



B. Traditional measures of phase coherence 

The classical approach to characterizing phase coherence 
of chaotic oscillations is based on the second-order structure 
function (variogram) of the detrended phase — (j){t) — 

2TTLUt, 



Di{.s) = (mt 



(10) 



Averaging this property over different realizations of the same 
process (or, as an alternative, over different time intervals cap- 
tured by the same trajectory - note that both options can be 
considered equivalent as long as the system under study can 
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be considered ergodic), one may approximately describe the 
dynamics of phase increments as a diffusion process Ii5ll6ll70l- 
l72l . In this case, one obtains: 

Dl{s) = B,s + Bo. (11) 

Comparing this with classical (stochastic) diffusion processes 
yields the phase diffusion coefficient D — Bi/2. We note that 
the proper estimation of this quantity from a single trajectory 
may be challenging, since the detection of a proper scaling 
window in which the above linear relationship holds may be 
a nontrivial task. This is particularly true for NPC dynamics, 
where the appropriate definition of the phase variable (f> is cru- 
cial. We note that the numerical values of D depend on which 
of the phase definitions from Sec. |III A| is used. 

As an alternative approach, in recent studies on the phe- 
nomenon of coherence resonance 1731 l74l . it has been sug- 
gested using the coherence factor 

CF - (T) /ar (12) 

(i.e., the coefficient of variation of Poincare return times Ti, 
with (r) and ctt denoting mean and standard deviation of T) 
as a measure of coherence of noise-induced oscillations. This 
approach can be directly transferred to the problem of distin- 
guishing PC and NPC deterministic-chaotic oscillations llTOl . 
given a properly selected Poincare section in phase space. 
However, we emphasize that in case of NPC chaotic oscilla- 
tions, the choice of such a Poincare section may be a difficult 
task itself. 



C. RP-based indicators of phase coherence 



On the other hand, since the empirical distributions of the lo- 
cal RN measures and log hy differ primarily with respect to 
their variance when comparing them for PC and NPC chaos 
(Fig. |5]l, we propose using the standard deviations ac and 
CTiog b as two further alternative measures for phase coherence. 
In addition, it may be helpful also considering higher-order 
statistics of the corresponding empirical distribution func- 
tions, e.g., their skewness 7c and 7iogb- 



IV. EXAMPLE I: BIFURCATION SCENARIO OF THE 
ROSSLER SYSTEM 

In order to systematically evaluate the performance of es- 
tablished as well as potential new RP-based indicators for 
phase coherence of chaotic oscillators, we study a part of the 
bifurcation scenario of the Rossler system (Eq. ([TJ), where the 
parameter a is systematically varied in the range [0.15, 0.3]. 
This parameter range comprises different kinds of dynamics, 
including periodic windows and PC as well as NPC chaotic 
oscillations. The transition between PC and NPC chaos oc- 
curs at flc ~ 0.2, which is in reasonable agreement with pre- 
vious studies using a slightly different parameter setting (e.g., 
Il68l ). Specifically, for a < Qc, the observed chaotic attractors 
are always PC, whereas they are NPC for a > Uc- In order 
to properly detect the location of periodic windows and sys- 
tematically exclude them when comparing the values of our 
measures for PC and NPC chaos, the largest Lyapunov expo- 
nents A1.2 of the system are additionally computed 19]. 



Since the proper estimation of the phase diffusion coeffi- 
cient D and coherence factor CF may be challenging, we 
will in the following use our results from Sec. |IIC| for defin- 
ing some novel indicators for phase coherence of chaotic os- 
cillations based on RPs. As we have already observed, the 
appearance of the RT distribution p(t) is different for PC and 
NPC chaos. Since mean RT (r) = (t) (e) and the corre- 
sponding standard deviation (Tt{£) do not provide sufficient 
results when considered separately, we suggest using the co- 
efficient of variation instead. This idea provides a straight- 
forward generalization of the coherence factor CF (based on 
the return times with respect to a fixed Poincare section) to 
a comparable measure based on the recurrence times to ar- 
bitrary e-neighborhoods of previously visited points in phase 
space. Consequently, we will refer to this measure as the gen- 
eralized coherence factor 



GCF = GCF{e) = 



(r) is) 



(13) 



Complementary to this approach, we also consider mea- 
sures characterizing the properties of the associated RNs. On 
the one hand, we suggest that some global network charac- 
teristics may be helpful for distinguishing between PC and 
NPC chaos, as they have already proven useful for discrimi- 
nating between complex periodic and chaotic orbits 11471 [50l . 



A. Traditional and recurrence times-based measures 

Figure [6] displays the variation of the Lyapunov exponents 
Ai_2, the phase diffusion coefficient D, the coherence index 
CI and the generalized coherence factor GCF when the pa- 
rameter a is changed. One clearly observes that the differ- 
ent measures are able to detect the transition between PC and 
NPC oscillations at about a — 0.21, but show different sig- 
natures in the presence of periodic windows. Specifically, 
the phase diffusion coefficient D takes values close to zero 
{D < 10^^) in both the periodic and PC chaotic windows, 
but gets much larger in the NPC chaotic regime (Fig. |6jJ). 
The latter observation coincides with a rather high variance 
for the NPC chaotic dynamics, which is mainly due to the 
subjectivity in choosing the scaling window for obtaining the 
linear regression parameters in Eq. (111. The coherence in- 
dex CI (Fig. [6p) is zero for a < 0.2, but strictly positive 
for higher values, including pronounced local maxima in the 
periodic windows (indicating that the periodic oscillations in 
these windows have no unique origin in the {x, ?/)-plane ei- 
ther). In contrast, the generalized coherence factor based on 
the recurrence time distributions takes very low values for 
NPC chaos, and higher ones for periodic and PC chaotic win- 
dows (Fig.|6]D). 
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FIG. 6: Behavior of different measures for phase coherence for the 
Rossler system in dependence on the parameter a (error bars indi- 
cate standard deviations obtained from 100 independent realizations 
of the system for each value of a): (A) Largest Lyapunov exponents 
Ai (solid line, A) and A2 (dashed line) calculated from the dynami- 
cal equations, indicating the location of periodic windows, (B) phase 
diffusion coefficient D, (C) coherence index CI, (D) generalized 
coherence factor GCF (RR — 0.03). Shaded areas indicate the 
presence of periodic windows evaluated by means of the largest Lya- 
punov exponents. 



B. Recurrence network measures 

In a similar way as described above, the values of local 
as well as global RN measures have been computed for re- 
alizations of the system for different values of a. Figure |7] 
shows the corresponding results. Regarding the global net- 
work characteristics, we find that the transitivity T has clearly 
higher values in the NPC regime in comparison to the phase- 
coherent chaos. In contrast, the assortativity coefficient TZ 
is clearly not capable of distinguishing both types of chaos, 
while a corresponding evaluation for C and £ requires more 
detailed statistical analysis (see below). Regarding the two lo- 
cal RN measures C„ and log 6„, standard deviation and skew- 
ness of both quantities show significantly higher values for 
NPC chaos than in the PC case, which is to be expected due 
to the more complex structure of the attractor in phase space. 
In general, the fluctuations of RN measures between different 
realizations obtained for the same value of a are much larger 
in the NPC regime than for PC chaos. For the periodic win- 
dows, T, C and C show pronounced maxima (which is con- 
sistent with previous findings B71 1501 ). whereas ac clearly 
displays local minima. In contrast, the signatures in TZ and 
the betweenness-based measure criog b are more complex. 



C. Discriminatory skills of RP-based phase coherence 
indicators 

In order to systematically compare the discriminatory skills 
of all proposed RP-based measures with respect to PC and 
NPC chaos, we divide the set of considered values of the con- 
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FIG. 7: Behavior of RN-based characteristics for the Rossler system 
in dependence on the parameter a (RR — 0.03, error bars indicate 
standard deviations obtained from 100 independent realizations of 
the system for each value of a): (A) global clustering coefficient C, 
(B) network transitivity T, (C) average path length C, (D) assortativ- 
ity coefficient TZ, and (E,F) standard deviation and (G,H) skewness 
of the local clustering coefficient and logarithmic betweenness cen- 
trality (ac, mogb, 7c and 7iog6, respectively). 



trol parameter a into three groups: one group Sq represent- 
ing the periodic windows (characterized by a maximum Lya- 
punov exponent Ai which does not significantly differ from 
zero within the numerical limits (i.e., Ai < A* = 0.02), and 
two groups Si and ^2 distinguished by whether or not the co- 
herence index CI (Eq. (j9])) does significantly differ from zero 
(i.e., CI{a) < Cr = 0.001 for PC chaos, and CI{a) > CI* 
for NPC chaos, respectively). Based on this initial discrimina- 
tion, we may statistically evaluate whether or not main statisti- 
cal characteristics of the distributions p{x\Si) of the different 
measures x obtained for both groups Si and 5*2 differ signifi- 
cantly. This problem can be solved by a one-way analysis of 
variance (ANOVA) |75|, with the factor being determined by 
two classes of values of CI. In order to evaluate whether the 
medians of some characteristic parameters in sets 5*1 and 52 
differ significantly (given the variances of the empiiically ob- 
served distribution functions), we perform a Mann- Whitney 
U-test |76t|77J, which can be considered as the equivalent of 
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0.84 (0.05) 


0.86 (0.04) 


0.2435 
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Clog 6 


0.56 (0.04) 


0.71 (0.09) 


1.18 X 10" 
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7c 


0.39 (0.52) 


-0.82 (0.62) 


1.33 X 10" 


11 




7log6 


-1.39 (0.65) 


-2.76 (0.48) 


8.47 X 10" 
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TABLE I: Mean values and standard deviations (in brackets) of the 
different measures for phase coherence for the considered realiza- 
tions of the Rossler system (averages over 100 independent realiza- 
tions for every value of a, fixed RR = 0.03) taken over all parameter 
values in the PC and NPC regimes, and P-values of the associated 
(7-test: generalized coherence factor GCF, global RN measures C, 
T, C and TZ, and standard deviation a. and skewness 7. of the distri- 
butions of the local RN measures Cv and log (from top to bottom). 
Symbols indicate the significance of the different parameters as dis- 
criminatory statistics ( — : insignificant, *: significant at 5% level, 
**: significant at 1% level, ***: significant at 0.1% level). 

an _F-test llTSll on the sets of rank numbers. 

The results of our corresponding analysis are summarized 
in Tab. |l] and confirm our qualitative statements. Specifically, 
we observe that standard deviation and skewness of the distri- 
butions of Cv and log 6„ allow a statistical discrimination of 
both chaotic regime with very high confidence. For the global 
RN measures, only network transitivity T performs compa- 
rably well. The average path length C also guarantees a re- 
liable discrimination, whereas global clustering coefficient C 
and assortativity coefficient TZ perform clearly worse. Finally, 
we find that the generalized coherence factor GCF based on 
the RT distributions in principle also allows distinguishing be- 
tween PC and NPC dynamics, however, on a much lower level 
of significance. 



D. Impact of the homoclinic point on RN measures 

A detailed inspection of the previously described findings 
reveals two interesting aspects: First, we observe that almost 
all RN-based measures show an overshooting close to the tran- 
sition between PC and NPC chaos (see Fig. |7]i. Further in- 
vestigations reveal that this effect does not result from a par- 
ticular choice of sampling or the finite length of the consid- 
ered realizations of the system (i.e., the presence of possibly 
transient behavior), but seems to be generic. Second, the be- 
havior of the network transitivity T seems to contradict re- 
cent general findings on the relationship between transitivity 
and effective attractor dimensions |50|: the higher the effec- 
tive dimension, the lower the RN transitivity. Specifically, the 
NPC regime has a higher dimension than PC chaos (this can 
be inferred from the higher maximum Lyapunov exponent in- 
dicating a higher Lyapunov dimension via the Kaplan- Yorke 
conjecture). Therefore, one has to expect that T takes higher 



values in the PC regime than for NPC chaos, whereas Fig.|7p 
displays the opposite behavior. (In a similar way, given the 
known fact that for a comparable value of e, periodic or- 
bits typically have a higher average path length than chaotic 
ones due to the formation of geometric "shortcuts" |45|, one 
would also expect C to be shorter in the NPC regime than for 
PC chaos, which is not consistent with the results in Tab. [l]) 
As we will argue below, these observations can, however, be 
explained in terms of the specific attractor geometry of the 
Rossler system, which is characterized by the considered RN 
properties. 

In order to understand the aforementioned overshooting as 
well as counter-intuitive behavior of RN measures, recall that 
the chaotic attractors of the Rossler system are characterized 
by the presence of a homoclinic point at the origin. In fact, 
the importance of the associated homoclinic orbit for the tran- 
sition between spiral-type (PC) and screw-type (NPC) chaotic 
oscillations has been widely recognized for the Rossler sys- 
tem as well as other chaotic oscillators with a similar tran- 
sition ||79H83]| . On the one hand, as the control parameter a 
increases within the PC chaotic regime, the attractor succes- 
sively grows and finally extends to the vicinity of the origin 
shortly before the transition to the funnel regime. On the other 
hand, the dynamics in the [x, ?;)-plane becomes very slow 
whenever a trajectory on the chaotic attractor gets close to 
the homoclinic point, before getting rapidly "ejected" out that 
plane following the direction of the associated unstable man- 
ifold. Thus, the growth of the chaotic attractor towards the 
origin has two consequences: First, the statistical properties 
of the distribution of ejection and re-injection "events" with 
respect to the (x, y)-plane changes markedly as a increases 
towards the transition point between PC and funnel regimes, 
which has a distinct effect on the overall recurrence properties 
of the system. This is reflected by the fact that the first return 
maps display one distinct differentiable extremum for spiral- 
type chaos, but several ones for screw -type chaos |83|. Sec- 
ond, due to the slow dynamics close to the homoclinic point, 
there is a high density of sampled points on a trajectory in the 
neighborhood of the origin, because the residence probability 
in this part of the phase space increases sharply shortly before 
the transition point. 

Within the framework of RNs, the accumulation effect 
around the origin becomes well expressed in terms of the dis- 
tribution of degree centrality ky = Tlj^v-^^.j, another im- 
portant local network measure. Specifically, while the mean 
degree (fc) — [N — 1)RR is constant when keeping the recur- 
rence rate (edge density) fixed (Fig. |8]\), the standard devia- 
tion increases strongly shortly before the transition between 
both chaotic regimes (Fig.[8j3), which implies the presence of 
many vertices with high degree, i.e., the existence of a phase 
space region with a high probability density of the attractor. 
As a result, the local network transitivity in this distinct region 
increases significantly: since the neighborhoods of many ver- 
tices (in our case those located close to the origin) are densely 
populated (high degree), they also show a high (local) cluster- 
ing coefficient (C^ < 1, cf. Fig.|3jj). This local behavior trans- 
lates into a higher global network transitivity T (Fig. |7]B) as 
well as a higher ac (Fig.|7]E). In a similar way, we can explain 
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FIG. 8: Mean values (k) (A) and standard deviations at (B) of the 
distribution of degree centrality fci, for tlie RNs obtained from 100 
independent realizations (error bars indicating ensemble means and 
standard deviations) of the Rossler system (TV = 10, 000). The de- 
sired recurrence rate RR — (k) /{N — 1) ~ 0.03 has been approxi- 
mated by selecting the threshold e based on a Monte Carlo sampling 
of inter-point distances from the trajectory in order to enhance com- 
putational efficiency. 
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TABLE II: As in Tab. |l] - results obtained with a fixed recurrence 
threshold e = 0.2776. 




FIG. 9: As in Fig. |7] for a fixed recurrence threshold e = 0.2776 
(corresponding to RR ~ 0.03 at a = 0.15). 



the overshooting of T and ctc close to the transition point, 
where the variance of the degree centrahty (and, hence, the 
density of points close to the origin) is the highest. 

Regarding the effect on the path-based measures £ and 
flog 6, we note that if we consider a fixed value of e (instead 
of a fixed RR) as a is changed, we find no overshooting close 
to the transition between PC and NPC chaos (see Fig. |9p,F) 
(in a similar way, the corresponding effect is clearly reduced 
for Clog 6 as well). Recalling the meaning of £ |451 . this ob- 
servation clearly indicates that the overall size of the attractor 
does not change markedly close to the transition point. More- 
over, £ now takes larger values for PC chaos than in the NPC 
regime (Tab.[Il)i, which reflects the increasing geometric com- 
plexity of the attractor In contrast to the path-based measures, 
the overshooting effect on the transitivity-based measures T 
and ac persists and becomes even enhanced for the global 
measures T and C, while it is reduced for ac and 7c . We em- 
phasize that with a fixed e, the recurrence rate RR becomes 
larger when increasing a close to the transition point due to 
the accumulation of vertices close to the origin, which could 
explain the aforementioned behavior. 



V. EXAMPLE II: BIFURCATION SCENARIO OF THE 
MACKEY-GLASS SYSTEM 

The scenario of a transition from spiral-type (PC) to screw- 
type (NPC) chaos is common to several nonlinear oscillators 
(for examples, see ll84l[85l ). However, there are further exam- 
ples for NPC chaos in other types of complex systems, espe- 
cially in time-delay systems. For illustrative purposes, in the 
following we reexamine a part of the the bifurcation scenario 
of the Mackey-Glass equation lIZTl 



xit) 



0.2x{t~T) 



0.1x{t), 



(14) 



a well-studied time-delay system, for r G [10, 20]. In this pa- 
rameter range, it is known that the system undergoes several 
transitions between periodic and NPC chaotic solutions ||28l 
(see Fig. 10 1. Note that unlike the Rossler system, the 



Mackey-Glass equation describes a time-delay system, i.e., an 
infinite-dimensional dynamical system. 

Figure [TT]\ shows the behavior of the maximum Lyapunov 
exponent with changing control parameter t. For r > 16, the 



10 



A 1.5 
1.2 
jL 0.9 

X 

0.6 
O.g 

C 1.5 
1.2 
jL 0.9 

X 

0.6 

0. 




0.6 



0.9 
x(t) 



1.2 



1.5 




1.3 



0.6 



0.9 
x(t) 



1.2 



1.5 



B 1.5 

1.2 
X 0.9 

X 

0.6 

°i 
D 1.5 

1.2 
X 0.9 

X 

0.6 
0.3 




0.6 



0.9 
x(t) 



1.2 



1.5 




0.6 



0.9 
x(t) 



1.2 



1.5 



FIG. 10: Phase portraits of the Mackey-Glass system l |14[ l for (A) 
r = 13, (B) T — 13.5 (after the period-doubling bifurcation), (C) 
T = 15.5, and (D) r = 17. 




10 12 14 16 18 20 TO 12 14 16 18 20 



FIG. 11: Behavior of different statistical characteristics for individ- 
ual realizations of the Mackey-Glass system in dependence on the 
parameter r: (A) Largest Lyapunov exponent Ai estimated from 
the variational equations of a discretized version of the system with 
10,000 variables representing {x{t), x{t - r/9999), . ..,x{t- r)), 
(B) phase diffusion coefficient D and (C) coherence index CI ob- 
tained via Hilbert transform of x{t). In addition, (D) shows the mean 
generalized coherence factor GCF obtained from 100 realizations 
for each value of r (RR — 0.03, embedding dimension 3 and delay 
r/2, i.e., X, = {x{ti),x{U - r/2),x{U - r))). 



Mackey-Glass system switches back and forth between peri- 
odic limit-cycle oscillations (Ai = 0) and chaotic solutions 
(Ai > 0). However, the phase diffusion coefficient D starts 
increasing from almost zero to non-zero (but still very small) 
values only at somewhat larger t (Fig. [TTj3), pointing to a 
gradual loss of phase coherence with rising control parame- 
ter. We note that this finding is distinctively different from 
those made for the Rossler system, where the system under- 
goes a rather sharp transition from PC to NPC chaos. The 
behavior of the coherence index CI (Fig. [TTp) based on the 




FIG. 12: Behavior of RN-based characteristics for the Mackey-Glass 
system in dependence on the parameter r (RR = 0.03, embedding 
parameters as in Fig. [TTp, error bars indicate standard deviations 
obtained from 100 independent realizations of the system for each 
value of a): (A) global clustering coefficient C, (B) network transi- 
tivity T, (C) average path length £, (D) assortativity coefficient TZ, 
and (E,F) standard deviation and (G,H) skewness of the local clus- 
tering coefficient and logarithmic betweenness centrality (ac, o-\og t, 
7c and 7iog b, respectively). 



standard Hilbert phase even shows a clear transition towards 
significantly positive values before the establishment of the 
first chaotic solution. This fact is clearly related to the spe- 
cific geometry of the attractor forming a small secondary loop 
structure after about r = 15 in the {x{t),x{t — T))-plane (see 
Fig. [Top). Finally, GCF (Fig. [TTp) shows a sudden drop at 
T > 13 (due to the presence of a period-doubling bifurca- 
tion 1 86 1 leading to marked changes in the RT distribution of 
the periodic solutions), followed by a clear downward trend 
for further increasing r. 

The above findings are further supported by the properties 
of RNs resulting from example trajectories obtained for dif- 
ferent values of r (Fig. \T2\ . As a first parameter interval 
of interest we consider t e [13, 14], which is characterized 
by Ai = 0, i.e., completely periodic dynamics. Here, all 
network measures show a marked transition indicating struc- 
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FIG. 13: Estimates of the marginal density p{x) of the Mackey-Glass 
system ^ for (A) r = 13, (B) r = 13.5, (C) r = 15.5, and (D) 
T — 17. One clearly observes the period-doubling bifurcation (B) 
and the emergence of a cusp point and, subsequently, the secondary 
loop structure (C) in terms of the present local maxima. 



tural changes of the underlying attractor corresponding to a 
period-doubling bifurcation. Specifically, C and T show a 
distinct drop from their values expected for periodic dynamics 
(C = T — 0.75 |50|), indicating the emergence of a structure 
of higher geometric complexity (cf. Fig. 10 \,B). A similar 



marked drop is shown by £, which is related to the emergence 
of geometric "shortcuts" after establishing the second loop of 
the periodic orbit. In contrast to C and T, this feature persists 
for higher t. In addition, TZ decreases abruptly at the period- 
doubling bifurcation. Regarding the local network properties, 
we find a sharp increase in the standard deviation, and a de- 
crease in the skewness of both clustering coefficient and log- 
betweenness distributions. We explain this observation by the 
fact that on the original single-loop limit cycle (with its rather 
homogeneous density), the local clustering coefficient does 
not vary much (crc ~ 0), whereas due to the emergence of 
the second major loop, there exists some "cross-over region" 
within which the neighborhood of state vectors has distinc- 
tively different shape and, hence, clustering properties in the 
recurrence network. It is interesting to note that at the same 
time, the associated skewness changes its sign as the period-2 
orbit successively develops (t < 14) before getting back to 
positive values for t > 14. 

A second interesting parameter interval with distinct 
changes of the attractor geometry is t G [14, 16], which 
still refers to the periodic regime of the Mackey-Glass sys- 
tem (Ai — 0). As shown in Fig.[TTp, at r w 15, the two-loop 
periodic orbit starts forming a cusp in the {x{t),x{t — t))- 
plane and, subsequently, an additional minor loop structure 
(Fig. 102), so that the associated Hilbert phase variable does 
not monotonously increase anymore. In parallel to this, for 
T > 14 both C and T increase beyond the expected values for 
a periodic orbit, although the maximum Lyapunov exponent 
clearly displays the presence of a limit cycle. This behavior 
results from an accumulation of probability on the trajectory 



close to the cusp (cf. Fig. [T3p). Since such an accumulation 
has a similar effect to a recurrence network as a fixed point 
(C„ = 1), the overall values of the transitivity-based measures 
increase. The same applies to C, where the presence of an ac- 
cumulation region leads to an overall reduction of the effective 
£ value to maintain the same recurrence rate RR. Similar con- 
siderations explain the increase of standard deviation and ab- 
solute value of skewness associated with the log-betweenness 
distribution. 

While so far only bifurcations between different periodic 
regimes have been discussed, for larger values of r we ob- 
serve that the transition from periodic to chaotic behavior is 
characterized by a sharp decrease in all four considered global 



recurrence network measures (Fig. 12 \-D), which is consis 



tent with the results obtained for the Rossler system |47| (see 
Fig. [?]) as well as other continuous -time dynamical systems. 
We note that unlike in other systems, both the periodic orbit 
prior to the transition point and the emerging chaotic solution 
are not phase-coherent (Fig.[TO]l. In this respect, the Mackey- 
Glass system does not allow studying geometric and dynamic 
differences between PC and NPC chaotic solutions, but serves 
as an illustrative example for the presence of noncoherent (pe- 
riodic and chaotic) oscillations and their impact on recurrence 
based characteristics. 



VI. CONCLUSIONS 

In summary, we have proposed a statistical framework for 
characterizing phase-coherent and noncoherent chaotic oscil- 
lations, which takes specific geometric information about the 
underlying attractor into account. For this purpose, we have 
utilized the recently developed concept of recurrence network 
(RN) analysis. Our results demonstrate that statistical mea- 
sures based on the recurrence properties of dynamical sys- 
tems do not only distinguish between periodic dynamics and 
chaos Il43ll47]| . or quasiperiodic dynamics and chaos ifTTHTOl , 
but also between different appearances of chaotic dynamics 
characterized by phase-coherent and noncoherent oscillations, 
respectively. In this spirit, RN analysis provides a widely ap- 
plicable tool for studying complex systems from a geomet- 
ric point of view, which supplements existing closely related 
techniques such as RQA and recurrence time statistics which 
characterize complementary properties directly related with 
the underlying dynamics. Specifically, besides studying sys- 
tems described by a finite set of ordinary differential equa- 
tions, it has been demonstrated that RN analysis is also ap- 
plicable for describing changes in the attractor geometry of 
time-delay systems such as the Mackey-Glass equation or the 
piecewise linear time-delay system studied in |87|. However, 
it is not yet possible to unequivocally distinguish between 
phase-coherent and noncoherent chaos exclusively based on 
individual characteristics of RNs such as network transitivity. 
The identification of a structural criterion for a corresponding 
discrimination will be a topic of our future work. 

For the Rossler system, we have studied the transition be- 
tween spiral- and screw-type chaos in some detail, which is 
common to several chaotic oscillators and leads to a change 
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from phase-coherent to noncoherent oscillations. The corre- 
sponding effects on the attractor geometry and, as a result, RN 
characteristics have been discussed in detail. As a particular 
result, we have shown that the recently given interpretation of 
the RN transitivity T, as a measure for the effective attrac- 
tor dimension, does not take statistical effects due to a very 
heterogeneous distribution of residence probability on the at- 
tractor into account, which has been largely overlooked in pre- 
vious research [50|. 

In general, we find that at least for the Rossler system sta- 
tistical characteristics based on the distributions of local RN 
measures allow an equal or even better discrimination be- 
tween phase-coherent and noncoherent chaos than some of the 
global network quantifiers. Moreover, both types of character- 
istics outperform the studied statistical characteristics based 
on the recurrence time distributions. We emphasize that RN 
measures probably behave so well because they explicitly 



characterize geometric attractor properties in phase space (i.e., 
no dynamic characteristics), which do strongly change at the 
transition between phase-coherent and noncoherent chaos. In 
this respect, they are particularly useful for obtaining a corre- 
sponding discrimination. 
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